High Magnetic Field ESR in the Haldane Spin Chains NENP and NINO 
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We present electron spin resonance experiments in tlie one-dimensional antiferromagnetic 5 = 1 
spin chains NENP and NINO in pulsed magnetic fields up to 50 T. The measured field dependence of 
the quantum energy gap for _B||6 is analyzed using the exact diagonalization method and the density 
matrix renormalization group method (DMRG). A staggered anisotropy term ( — l)'d(S'f Sf + 3^3^) 
was considered for the first time in addition to a staggered field term { — lySfB^^. We show that 
the spin dynamics in high magnetic fields strongly depends on the orthorhombic anisotropy E. 
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According to the Haldane conjecture W the ground- 
state properties and the excitation spectrum of a one- 
dimensional Heisenberg antiferromagnet are determined 
by the spin value 5* in a fundamental way. For inte- 
ger spins the system exhibits a disordered nonmagnetic 
groundstate which is separated from the first excited 
triplet states by an energy gap, whereas half-integer spin 
systems have no gap. The existence of the Haldane gap 
has been experimentally confirmed by inelastic neutron 
scattering (INS) experiments [^1 and magnetization mea- 
surements ^. 

Electron spin resonance (ESR) experiments have 
proven to be an excellent method to study the mag- 
netic field dependence of the excitation spectrum in 
Haldane systems Q - [|o|. In the prototypical 
one-dimensional spin chains Ni(C2H8N2)2N02(C104) 
(NENP) and Ni(C3HioN2)2N02(C104) (NINO) transi- 
tions between excited states have been observed at the 
Brillouinzone boundary at q — n/b nM- In addition 
transitions from the singlet groundstate at g = to the 
excited triplet at q = n/b have been found 0- pffl, which 
are normally forbidden by selection rules. These ESR 
experiments could be explained using the concept of a 
transverse staggered magnetic field [nill2| which mixes 
the eigenstates at q = and q = ir/b. Such a staggered 
field could be induced by an external magnetic field due 
to the alternating tilting of the local anisotropy axes of 
the Ni^+ ions and was first observed in NMR experiments 
|l3| . ESR experiments are so far the only method for a 
direct observation of the energy gap in Haldane systems 
well above the critical field B^ (8.5 T in NENP and w 8 T 
in NINO). 

NENP and NINO both crystallize in the orthorhom- 
bic system and show similar lattice constants [Q - f[l6| . 



NENP can be described in the Pnma space group Iff:] 
whereas NINO belongs to the Pbn2i space group ||lj,|l6[. 
The Ni^+ ions are covalently linked via N02-complexes 
which are responsible for the antiferromagnetic superex- 
change interaction. Both compounds are good realiza- 
tions of one-dimensional chains since the interchain ex- 
change J' is about a factor of 10'* smaller than the in- 



trachain exchange J |l4| due to the separation by the 
CIO4 ions. The local surrounding of the Ni^+ ions is a 
distorted octahedron with a basal plane built out of four 
nitrogen atoms from the diamine groups. The axial po- 
sitions are occupied by an oxygen and a nitrogen atom 
fronr the nitrite bridges. This local surrounding gives 
rise to a planar anisotropy D and to an orthorhombic 
anisotropy E as well as to an anisotropic g tensor. 

A significant characteristic in NENP and NINO is 
the alternating tilting of the local surrounding of the 
Ni^"*" ions along the chain. As we will show in this pa- 
per, it is important to know the direction of this tilting 
with respect to the easy direction of the orthorhombic 
anisotropy E. For NENP it is clear from the structure 
||l5| and NMR experiments ||l|l that the tihing of the 
local environment is along the c direction p7[ while the 
easy direction of the orthorhombic anisotropy is along a 
(see e.g. Q). In NINO the easy direction in the basal 
octahedron plane is along c [Q while the structure [l^ ] 
gives rise to the assumption that the local environment is 
tilted towards the a-axis. However it is not clear how to 
deduce a definite angle of tilting of the g or the D tensors 
out of the structure. 

The ESR experiments were performed in the high mag- 
netic field laboratory at the University of Frankfurt. Dif- 
ferent single crystals of NENP and NINO have been 
mounted in a transmission type sample holder in Voigt 
geometry (propagation vector k perpendicular to the ex- 
ternal magnetic field B) with the chain direction b par- 
allel to B. We used Gunn and IMPATT diodes together 
with frequency multipliers to produce mm- wave radiation 
in the frequency range 55-200 GHz. The radiation was 
unpolarized when passing the sample and was detected 
with a fast InSb hot electron bolometer mounted in a 
separated cryostat. High magnetic fields were produced 
using an AI2O3 strengthened copper coil with internal 
reinforcement, manufactured in the National High Mag- 
netic Field Laboratory in Tallahassee, Florida. Fields 
up to 50 T with 8 ms rise time could be produced by 
discharging capacitors with a stored energy of 0.4 MJ. 
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FIG. 1. Transmission spectra in NENP and NINO at 
T = 1.5 K. The arrows mark excitations from the groundstate 
below and above the critical field in the 134 GHz spectrum. 
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FIG. 2. Frequency field dependence of the groundstate ex- 
citations in NENP and NINO (full and open symbols respec- 
tively) . 



While the coil is operating at 77 K, the sample was 
cooled down to 1.5 K in a pumped ^He cryostat. The 
magnetic field was measured with a pick up coil and cal- 
ibrated using a DPPH reference sample. Details of the 
high field setup can be found in |19 

Figure m shows some typical mm-wave transmission 
spectra at r=1.5 K for NENP and NINO in high mag- 
netic fields. A drastic increase of the resonance inten- 
sity can be observed with increasing field above the crit- 
ical field Be ~ 8 T (see corresponding arrows in fig. H). 
This increase of intensity was explained before in the 
framework of the staggered field concept 0,p2l. Also 
the linewidth increases remarkably in both substances: 
while the typical linewidth for B < Be in NENP and 
NINO is about 0.1 T and 0.05 T respectively, it increases 
for B > Be about a factor of 20 in NENP and a factor 
of 10 in NINO. No decrease of the integrated intensity 
could be observed with decreasing temperatures, indi- 
cating that the observed transitions are excitations from 
the groundstate y|. The resonance frequencies therefore 
give a direct measure of the magnetic field dependence 
of the energy gap ||ll|,|l||. 

The frequency field dependence for all observed reso- 
nances is shown in fig. g together with low field data ob- 
tained by other groups. For B < Be the gap in both sys- 
tems decreases with a slope corresponding to the g factor 
(gt = 2.15 in NENP, gt = 2.17 in NINO). The gap does 
not close but remains finite at the critical field Be and 
increases again for B > Be- For B > Be the curvature of 
the A{B) curve is much stronger in NENP than in NINO. 
The smaller slope in the high field region for NENP seems 
to be one explanation for the larger linewidth in this com- 
pound. 

The difference in the high field spin dynamics between 
NENP and NINO was thought to be mainly caused by 
a larger staggered inclination of the Ni^+ surrounding in 
NENP |^,|l^ (although these assumptions are based on 
calculations with an isotropic Hamiltonian and unrealis- 
tically small exchange constants J). Thus the ratio of 



the gap energies and the ratio of the slope of the A(B) 
curves are expected to be approximately constant be- 
tween Be and B = 4 J 1 11 1 . In contrast the actual field de- 
pendence of the gap energies shows a different behaviour 
with a crossing point of both A(B) curves at i? = 1.8J 

(fig i- 

To analyze our data for the B\\b configuration more 
realistically we use the following Hamiltonian 



H = Y,{JS^S^+l + D{Slf + ^[(5f )2 

i 

- g.lM^StB + {-md{S^S! + StS^) + 



cg,^,BStB]}. (1) 



The indices x, y and z correspond to the crystallographic 
axes c, a and b respectively. J is the antiferromagnetic 
intrachain exchange, D and E are the planar and the 
orthorhombic anisotropics. The fourth and the follow- 
ing terms contain the Zeeman coupling and the staggered 
terms due to an alternating tilting of the anisotropy- and 
the g-tensor within the x-z plane, c is the staggered field 
constant c = Bst/B which is given by the tilting angle 9g 
and the difference of gx — gz'- 



gx- gz l 

gz 2 



tan20n 



(2) 



d is the staggered anisotropy constant which is given by 
Od and the difference oi D ~ E: 



d= {D- E)-tSin20D 



(3) 



The parameters D, E, g^ and gz are the observable val- 
ues in the crystallographic system. They can be derived 
by averaging over the tilted anisotropy- and g-tensors of 
neighboring Ni^+ sites. For a more detailed derivation of 
(^ and (||) we refer to ref. Q . 

We like to point out that although the Hamiltonian dl^) 
contains numerous parameters, only the angles 9g and Ou 
can be considered as " free" parameters within the fitting 
process for our high field data. J, D and E are fixed by 



0.30 
0.25 
0.20 

,0.15 

3 

0.10 
0.05 



0.00 







' ' ■ 










\ 












<^' 


"\\ 






^^ 


'^' 




^^^ : 




1 


5^-' 


"^ 


O/ 


• 


NENP : 




¥ 


^ 






o 


NINO 


■ 


1 




1 , 1 


1 


1 


1 , 1 , ■ 



0.0 0.2 0.4 0.6 0.8 



1.0 
B/J 



1.2 1.4 1.6 1.8 2.0 



FIG. 3. Magnetic field dependence of the gap in NENP and 
NINO in units of J. Symbols are experimental data, lines are 
extrapolated numerical results based on the Lanczos method 
(for B/J < 0.4) and the DMRG method (for B/J > 0.4). 
The staggered field and anisotropy parameters are c=-0.026, 
d/J=0.06 and c=0.0105, d/J=0.06 for NENP and NINO re- 
spectively (solid lines). The dashed line is the best fit for 
NENP with d/J = using c = 0.017. 

the zero field size of the gaps Ai, A2 and A3 between 
the groundstate and the excited triplet states a,t q = ir/b. 
J is given by the well known ratio 0.4U = Eq — (Ai + 
A2 + A3)/3 Q while D and E can be deduced from 
the splitting of the triplet A3 - (Ai + A2)/2 ^ and 
A2 — Ai respectively. The gap values are taken from 
ESR experiments |9[| which are in good agreement with 
INS results |2|. The g- factors for the different directions 
are taken from zero field susceptibility data ( [|5| for 
NENP and |l|l for NINO). 

Accordingly we obtain the following parameters 
J=43 K, D/J=0.17, E/J=0.01, g^=2.21, and 5^=2.15 
for NENP and J=39 K, D/J=0.2, E/J^-0.027, 5^=2.23, 
and 5^=2.17 for NINO. As noted before, in NINO the 
easy direction of the orthorhombic anisotropy is along 
the c- (x-) axis leading to a negative sign of E. Since the 
alternating tilting in NINO is supposed to be in the y-z 
plane the Sf operators and gx in the staggered term of 
(^ have to be replaced by Sf operators and gy respec- 
tively. 

The magnetic field dependence of the gap was calcu- 
lated both with exact diagonalization for A^=16 spins 
using the Lanczos algorithm and with the density matrix 
renormalization group (DMRG) method. The Lanczos 
results for short chains can be extrapolated in a reliable 
way for B < Be- For B > Be however the gap shows 
large finite size effects which manifest themselves in os- 
cillations. This is in particular the case for small values of 
the staggered fields, when both the ground state energy 
and the first excited state energy follow almost polygonal 
lines as a function of B. Therefore for B > Be we use the 
DMRG method which allows to treat large systems with 
high accuracy by keeping only the most important eigen- 
states of a blocked density matrix constructed according 




FIG. 4. Magnetic field dependence of the high field energy 
gap for different parameters. Full circles (•) are calculation 
for c = 0.0105, d/J = 0, D/J = 0.2 and E/J = 0.027. Open 
symbols are calculations which differ from (•) in only one 
parameter ((o) c = 0.012, (v) d/J = 0.06, (o) E/J = 0). 
The inset demonstrates the linear dependency of the gap on 
the orthorhombic anisotropy E at B/J = 2 for c = 0.0105, 
d/J = 0.06 and D = 0.2. All results are extrapolations of 
DMRG calculations. 

to ref. pj. We found best convergence using DMRG 
for N — 40 . . . 120 spins with open boundaries but also 
employed periodic boundary conditions in order to re- 
move the S = 1 artefact state occurring for open bound- 
aries in the low field range. In the high field region the 
method converges very well with as little as 30 states. In 
the region near the critical field more states and larger 
systems arc needed. 

Figure shows the extrapolated results of our calcula- 
tions together with the experimental data in units of J. 
For NINO we find c = 0.0105 and d/J = 0.06 assuming a 
tilting angle ofOg = 60^ 17.5°. For NENP we obtained 
our best fit for c = —0.026 and d/J — 0.06 based on the 



tilting angles 6g 



-31° and Od = 18.5° 



which differ 



in size and direction. 

This somewhat strange result implies that the the g- 
and the D- tensors in NENP depend on different neigh- 
boring atoms of the Ni^+ ions, e.g. one tensor is mainly 
determined by the position of the nitrogen atoms of the 
ethylenediamine groups whereas the other tensor is prin- 
cipally given by the position of the nitrite bridges. A 
more detailed analysis is needed at this point to clarify 
the dependence of the g- and D-tensors on the Ni^+ en- 
vironment. Nevertheless neglecting the staggered d term 
leads to a rather poor fit for the new high field data of 
NENP (dashed line in fig. m showing that this d term 
is essential to describe the gap as a function of the field. 
Besides one has to keep in mind that the interchain ex- 
change J' was not taken into account in our calculations 
which possibly could affect the high field spin dynamics. 

It is of special interest to study the effect of the differ- 
ent terms in (H^) on the energy gap when varying several 
parameters. In a first step we checked that finite values 



of c and d have practically no effect for B < Be, thus 
confirming the values taken for J, D, and E. Next we 
studied the effect of c and d in the high field range. For 
B > Be the system is expected to be gapless without any 
staggered term in (||). The staggered terms lift the de- 
generacy between the magnetic groundstate and the first 
exited state in the high field phase yielding a gap which 
increases with B. Fig. H shows the gap for different c and 
d values as a function of field. The effect of the d term 
on the gap is quite small at the critical field Be ~ 0.4J 
but increases with B comparable with the increase of the 
magnetization J2 < ^z > reflecting the effect of the S^ 
operators in (p. 

As is well known the anisotropics D and E reduce the 
gap for B < Be [£lj , |§2| . While the planar anisotropy D 
only has a very weak effect for B > Be it turned out that 
the orthorhombic anisotropy E considerably changes the 
size of the high field gap induced by the staggered c and 
d terms (fig. |j with inset). The gap is linearly reduced 
for a positive E in (JlJ) whereas it is linearly enlarged 
for E < (see inset of fig. ^). The staggered c and d 
terms induce a staggered magnetization J2i~^T < ^f > 
along the x direction in NENP ^^ (along the y direction 
in NINO), which is the hard axis of the orthorhombic 
anisotropy respectively. Thus a positive E increases the 
energy of the groundstate and the first excited states. 
This increase is somewhat stronger for the groundstate 
which shows a higher staggered magnetization than the 
first excited states yielding a reduction of the gap. 

The orthorhombic anisotropy term is essential to de- 
scribe the differences in the high field spin dynamics of 
NENP and NINO. As shown in fig. | the field depen- 
dence of the gap has a clear curvature for _E = while it 
becomes more linear with increasing E corresponding to 
the experimental results in both substances. 

In the present work we have shown unprecedented 
data for the quantum energy gap in the Haldane sys- 
tems NENP and NINO measured by ESR experiments in 
pulsed magnetic fields up to 50 T. We presented numeri- 
cal calculations based on a direct diagonalization and the 
DMRG method considering the planar and orthorhombic 
anisotropics D and E as well as the staggered anisotropy 
and the transverse staggered field. The high field data 
could be fitted using realistic values for the exchange in- 
teraction and the anisotropics. In addition to a stag- 
gered field term the staggered anisotropy term is essen- 
tial to describe the curvature of the gap as a function of 
field. Finally the orthorhombic anisotropy E turned out 
to have an unexpected strong effect on the energy gap 
for B > Be. 

Future experiments with high magnetic fields perpen- 
dicular to the chain axis could provide us with more infor- 
mation especially for the case where no staggered field is 
induced. In addition experimental and theoretical inves- 
tigations on the field dependence of higher excited states 
should extend the knowledge about the high field phase 



in these Haldane systems. 
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